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Experiments on tetragonal Sr2Ru04 and Celrlns indicate the presence of superconductivity with 
a multi-component superconducting order parameter. Such an order parameter should exhibit an 
in-plane anisotropy in the upper critical field near the superconducting transition temperature that 
\ does not occur for single-component superconductors. Here, this anisotropy is determined from 

microscopic calculations for arbitrary gap functions. It is shown that this anisotropy is generally 
' not small and, in some cases, independent of impurity scattering. Furthermore, this anisotropy 

is calculated for many detailed microscopic models of Sr2Ru04. For these models this anisotropy 
is found to be large, which is in sharp contrast to the small anisotropy observed experimentally. 
OO ' However, an accidental cancellation of the anisotropy for gaps on different Fermi surface sheets can 

lead to a result that is consistent with experiment. 
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^ ■ One of the central issues in the field of unconventional superconductivity is the identification of the order parameter 
O ' symmetry. This issue can usually only be [addressed by determinations of the phase of the order parameter (such 
^ . as through Josephson junction experimentstJ) . However, in the case that the Cooper pair wave functions have more 
' than one complex degree of freedom, other experimental signatures can reveal the pairing symmetry. In this paper 
one of these signatures, the upper critical field anisotropy for fields in the basal plane, is explored for tetragonal 
^ • superconductors. This anisotropy was established by Gor'kov on phenomenological groundso. Such an anisotropy, for 
1^ \ which dHc2/dT\T=T^ is not equal for the magnetic field apphed along the (1,0, 0) (along the crystallographic a axis) 
' and the (1,1,0) directions, cannot occur for superconducing order parameters that have only one complex degree 
^ ■ of freedom (for example conventional s-wave superconductors or (ia;2_j,2-wave superconductors). Consequently, this 
provides a clear test of the pairing symmetry. For tetragonal materials, there exist two possible pairing states that 
exhibit this property. One is the spin triplet i?„ representation which can be described a gap function of the form 
(note that this is not the most general form) d(k) — zlrjifxCk) +ri2fy{'k)], where in the simplest case (for a cylindrical 
Q ' Fermi surface,) /a;(k) = kx and /y(k) = ky. Such a gap function is believed to describe the pairing state in Sr2Ru04 
I— ' I (see below )l3'Q. The other pairing state is the spin-singlet Eg representation which can be described by a gap function 
of the form ^'(k) = r]ifxz(k.) + il2fyzO^), where fxzO^) = kxkz and /j,z(k) = kykz in the simplest case. In this article, 
^ the anisotropy is determined within microscopic theories with arbitrary gap functions. It is shown that the anisotropy 
. is generally not small and consequently easily observable. Furthermore, it is shown that for the Eg representation 
' and for some theories for the Eu representation, this anisotropy is impurity independent. Application of these results 
lO . to Sr2Ru04 and Celrlns are discussed and detailed calculations for Sr2Ru04 are also presented. Prior to discussing 
the anisotropy, an overview of the superconductivity in Sr2Ru04 and Celrlns is given. 

The oxide SX2RUO4 has a structure similar to high Tc materials and was observed to be superconducting by Maeno 
\^ et al. in 1994B. It has been established that this superconductor is not a conventional s-wave superconductor: NQR 
I measurements show no indication of a HebeL-Slichter peakQ in 1 /TiT, and Tc is strongly suppressed below the maximum 
. value of 1.5 K by non-magnetic impuritiea3. More recent experiments indicate &iy,Eu odd parity gap function of the 
form d(k) = z[riifx{}i) + 772/jf(k)]. The Knight shift measurements of Ishida et alu reveal that the spin susceptibility 
^ is unchanged upon|^ntcring the superconducting state; this is consistent with p-wave superconductivity (as predicted 
'■^ by Rice and Sigristlj and Baskaranllj). Furthermore, these measurements were conducted with the applied field in the 
^ basaLplane. Since the orientation of the gap function (that is d) is orthogonal to the spin projection of the Cooper 
O pairM, these measurements are consistent with the gap function aligajed along the z direction. The /zSR experiments 
of Luke et. al. have revealed spontaneous fields P| Meissner stateEj. This indicates that the superconducting order 
^ : parameter must have more than one componentOll3. This leads naturally to the conclusion that the superconducting 
gap function has the £orm d(k) = z{fxiify). This conclusioa is consistent with measurements of the field distribution 
of the vortex latticeEJ, Josephson tunnelling experimentstj, and point-contact measurementsEj. One immediate 
■ consequence of such a gap function is that the resulting low energy excitation spectrupuis gapless in the cleanJimit. 
However, recent experiments indicate thisnis not the case. Specific heat measurementsEj, NQR measurementstZI, and 
cavity penetration depth measurementsES all indicate the presence of nodes in the superconducting gap function. 
This has prompted a recent series of proposals that the gap function has the fO|rBi,d(k) = z[r7i/a;(k) -I- ry2/n(W] where 
/a;(ko) = /i/(ko) = for a set of points kg that form a line on the Fermi surfaceEj^Eil or are very anisotropic&Ej. These 
gap functions are all formally of the same E^ symmetry. The fact that these theories all have the same symmetry 
implies that there exist experimental signatures that will be common to all of them. In particular, the uppieij critical 
field anisotropy discussed above should exist. Recently, such an anisotropy has been observed by Mao et alBB. It was 
determined that 1 — i?c2"°V^c2^"' = 0-03 at a temperature of 0.35K (« Tc/5), 1 — ^^c2°"V-^c2^°'' observed to 
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change sign as T —f Tc, and appears to remain small for T near Tc- In this paper, the anisotropy near Tc will be 
calculated for a variety of models. 

Very recently, it has been reported that . /.(SR experiments of Heffner et. al. have revealed spontaneous fields 



in the Meissner state of tetragonal CelrlnsEU'ES. Just as in the case iQp-Sr2Ru04 above, this indicates that the 
superconducting order parameter must have more than one componenttillla. This leads naturally to the conclusion 
that the superconducting gap function is either of symmetry of of Eg symmetry. Consequently, both pairing 
symmetries are considered in this article. i— . 
The free energy for the Eu and Eg representations of 1)4/1 is given bylHl 

•12 I o 1^14 10 \ R l„ ^* ^ ^*^2 , 



/ = - a|^|^ +/3i|77l72 + /32(^ir?2* -^72^i*)72 



+ ^^2{\Dym? + \D^M'') (1) 

+ K3[{D,r]i)iDyr,2)* + h.c] 
+ i^A[{Dym){D,r)2y + h.c] + /iV(8^). 

where Dj = Vj — iff^jj h = V x A, and A is the vector potential. The stable homogeneous solutions are easily 
determined0. There are three phases: (a) rf = (l,i)/%/2 {(32 > and /32 > P3/2), (b) fj = (1,0) {(33 > and 
P2 < /33/2), and (c) fj = (1, 1)/V2 {(33 < and P2 < 0). The weak-coupling approximation predicts that phase (a) is 
stable since this phase minimizes the number of nodes in the order parameter [note that in exceptional circumstances 
phase (a) may be degenerate with phase (b) and (c) in the weak-coupling limit, for example in the theory of Ref.E3 
for = 1]. 

an external magnetic field applied in the basal plane, the upper critical field can be determined and is found to 

. , ahc , , 

Hc2{0) = = (2) 



e^/2^y i^2 + Ki~ \l (ki - K.2Y cos'^{29) + (K3 + sin2(2^ 



where 9 is the angle between the magnetic field and the crystallographic a axis. From the above expression, the 
ratio of Hc2 for the field along the (1,0,0) direction to that for the field along the (1,1,0) direction is given by 
^U^i.o)^jjU.o,o) _ ^ j.^^ _j_ ^2 — -|- K4\)/2min{K2, ki) where min{K2, ki) means the minimum |Of ki and K2. Note 
that Sigrist has also determined an additional anisotropy that occurs due to spin-orbit couplingEJ. This additional 
anisotropy introduces a correction to Eq. |^ that vanishes as T — > T^, so it is not considered here. 

Here, this anisotropy is determined in the weak-coupling approximation. This approximation is reasonable for 
Sr2Ru04 since Tc/Tp « 10"'', however it is not clear whether or not this approximation is valid for Celrln5. Using 
a gap function of the form d(k) = z[r]ifa:{\i) + 'i]2fy{\i)] for the Eu representation and introducing isotropic impurity 
scattering within a T-matrix approximation yields (the details will be given in a later publication) 

— 2 l^n>0 (w„+r)^ [Jx^y/ 

'^3 — 2 Z^n>0 (w„+r)3 ^ 

[{fxfyVxVy) + ^^^ifxVx)^] 

7rJV(0) 1 Iff \ 

^^4 — 2 2^n>0 (w„+r)^ KJxJyVxVy) 

1 

>o (w„+r)a 



7rAf(0)v^ 1 /f2„,2 



where aj„ = {2n + l)7rT, A^(0) is the density of states at the Fermi surface, F = I/tat is half of the scattering rate, Vi 
are the components of the Fermi velocity, and the brackets () denote an average over the Fermi surface. It has also 
been assumed that (/^ -I- fy) — 1. The terms proportional to F in ki and K3 are vertex corrections (which have the 
same form as those that arise in the calculation of current-current correlation functions). In the clean limit, this gives 



l {flvl) + {flvl)~mxfyVxVy)\ 

Hg^''''> V 2rmn{{f^vl),{f^vl)) ' ^' 
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For K2 < Ki and K3 + K4 > 0, this ratio is independent of the impurity concentration. The requirement that K2 < ki 
impUes for the field along (1,0,0) that (?yi,r?2) = ??(r)(l,0) [for K2 > ki (?7i,r?2) = ??(r)(0,l)] while the requirement 
that K3 + K4 > implies that for the field along the (1, 1, 0) direction that (r/i, 772) — ?7(r)(l, 1) [if K3 + K4 < then 
(??i7?72) = '7(r)(l, —I)]- Note that as the magnetic-^eld is decreased, there will exist a second order phase transition 
due to a change in the order parameter structureEd. For example, for the field along (1,0,0) and for K2 < ni there 
will be a second order transition with decreasing field where the component (771,772) oc (0,1) ^ 0. Such a transition 
reduces the number of nodes in the gap structure which is why it occurs. 
For the Eg representation using \l/(k) = {rjif^z + 112 fyz) gives 

7tN{0) 1 /f2 „.2\ 

'^l ~ 2 l^n>0 {uj„+r)'-'\Jxz''x/ 

7tN{0) 1 /f2 ,.2\ 

^2— 2 Z^n>0 (w„+r)a Uxz^/ 



TrAf(O) 1 



'^3 — 2 Z^n>0 
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[(/xz/yz VxVy) (4) 

ttN(0) 1 / f f \ 

'*4 — 2 2^n>0 (w„+r);* \JxzJyzVxVy) 

7rjV(0) 1 /f2 ,2\ 

'^S — 2 2^n>0 (w„+r)a KJxz^zI 

The main difference with the -En representation is that no vertex corrections appear. For an arbitrary impurity 
concentration this gives 



Hg^°^'^ V 2rmn{{P,vl),{PzVl)) ' ^> 

This ration is independent of the impurity concentration. These results (and the more detailed calculations below) 
indicate that the anisotropy is not generally small and should therefore be easily detected. It would be of interest to 
look for this anisotropy in Celrlns. 

It is informative to determine the anisotropy for existing microscopic theories in the clean limit. Since no such 
theories for Celrlns exist to date, the remainder of the analysis concentrates on Sr2Ru04. To do this the electronic 
dispersion must be given-asjiiust the form of the gap function. Initially, consider theories that for simplicity use a 
cylindrical Fermi surfaceEJod'Ej. For a gap function with nodes in the basal plane {fx,fy) — [(fc^ ~ ^y)^x,—(,k'^ — 

fc^)fcy]EB0, nj^l^^^ / — 0.58. For an anisotropic but fully gapped gap function: {fx,fy) = [sin(fc^), sin(fc2^)] 

with kp = 0.97tEj, H^l'^'^^ / H^l'^'^^ = 0.53. These estimates indicate that the anisotropy should be easily observable 
if the gap contains nodes in the basal plane or ia-Kecji anisotropic (these possibilities might be expected due to 
recent experiments that indicate nodes in the gapll3^ll3) . This in contrast to the experimental results of Mao et ali 
A more careful investigation below indicates that the anisotropy can accidentally be hidden, even for theories with 
nodes in the plane. 

To consider more realistic models, a tight binding approach is used to describe the electronic dispersion near the 
Fermi surface. Local-density approximation band-structure calculations r av a al that the density of states near the 
Fermi surface is due mainly to the four Ru 4d electrons in the t2g orbitalsBc3. There is a strong hybridization of 
these orbitals with the O 2p orbitals giving rise to antibonding tt* bands. The resulting bands have three quasi-2D 
Fermi surface sheets labelled a, (5 and 7 (see Fig. 1). To model these Fermi surface sheets, the following tight-binding 
dispersions are used: 

= e° — 2i.^(cos kx + cos ky) — At^ cos kx cos ky 

Co 



= e°_^ - ta.picoskx + COS ky) ± ^Jt\p{coskx - cosfcy)^ -|- IGi^^^sin^ kx sin^ ky (6) 

In following tight binding values are used: [e^,t^,t^) = (—0.4,0.4,0.12) for the 7 sheet and the values 
(Eq /3: ia,/3, ia,/?) — (—0.3,0.25,0.0375) for the {a,/3} sheets0. Table 1 shows the anisotropy ratio for gap functions 
found in various theories for Sr2Ru04. 

Note that for the a and (3 bands the gap functions with Vx replaced by sin kx were not included (as they were for 
the 7 band) . This is due to quasi one-dimensional nature of the a and (3 sheets. In particular, Vx.a oc sinfc2:[l -t- 
sgn{coskx — cosfcj^)]/2 ^ sinkx in the one-dimensional limit. It is Vx, not sinfe^;, which would be the gap functioa 
for a one-dimensional nearest neighbor pairing interaction. This is-supported by the calculations of Kuroki et al.ci 
where the {fx, fy) = {vx,Vy) basis functions on the a sheet of Ref.Eil gives a good approximation to the crib-shaped 
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gap function found in this reference [the {fx,fy) = {sin k^, sin ky) basis functions gives a quaHtatively different gap 
function, with gap maxima along the (1,1,0) momentum direction]. A point of interest is that 11^2^ /Hl^^ < 1 for 
the all the 7 sheet gap functions and Hl2^ /HH^ > 1 for all the a and /3 sheet gap functions. Consequently, the 
upper critical anisotropy can cancel between these sheets of the Fermi surface. For example, using the gap function 
fx = Vx and taking the magnitude of this gap function on the a and (3 surfaces to be the same (ca = cp) the gap 
magnitude ratio c^/ca — 6.7 will lead to zero upper critical field anisotropy (where the experimentally measured 
ratios N^{Q) : N(3{0) : iVa(O) = 1.334 : 0.914 : 0.216 and vj.^^ : vj.^^ : ^ = 2.6 : 6.8 : 7.8 were usedB). The same 
can be done for gap functions with line nodes: for fx = VySinkxSinky^ c-^/ca — 1.7 and for fx — {coskx — cosky)vx, 
Cy/ca — 55. These estimates imply that for such an accidental cancellation to occur the maximum gap on the 7 Fermi 
surface sheet must be larger than that on the {a, /?} sheets. Note that in principle all the gap functions listed in Table 
1 can appear simultaneously in the pairing state of Sr2Ru04 since they are all of the same symmetry. 

The in-plane anisotropy in the upper critical field near Tc has been calculated for microscopic models of the Eg 
and Eu pairing symmetries of the tetragonal point group. It has been shown that this anisotropy should be easily 
observable and may be present in Celrln5. Also it is shown that for Eg gap functions and for some E^ gap functions, 
this anisotropy is independent of impurity scattering. For Sr2Ru04, it is shown that an accidental cancellation of this 
anisotropy can occur if the 7 sheet of the Fermi surface is largely responsible for the superconductivity. 

I wish to thank V. Barzykin, L.P. Gor'kov, R.B. Joynt, Y. Maeno, M. Sigrist, and L. Taillefcr for informative 
discussions. 
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FIG. 1. Fermi surface of Sr2Ru04 using the tight binding parameters in this paper {kx and ky are in units tt/o where a is the 
lattice spacing). Moving out from the k = (0,0) point, first the /3 sheet, then the 7 sheet, and finally the a sheet are crossed. 
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TABLE I. Upper critical field anisotropy for various gap functions for Sr2Ru04. Note that the normalization of the gap 
functions has not been included here. The column called Impurity Dependent states whether or not Hl2° /HI2'' depends upon 
impurity scattering. 
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